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We investigate the effect of the trapping potential on the quantum phases of strongly correlated 
ultracold bosons in one-dimensional periodic and quasiperiodic optical lattices. By means of a 
decoupling meanfield approach, we characterize the ground state of the system and its behavior 
under variation of the harmonic trapping, as a function of the total number of atoms. For a small 
atom number the system shows an incompressible Mott-insulating phase, as the size of the cloud 
remains unaffected when the trapping potential is varied. When the quasiperiodic potential is 
added the system develops a metastable-disordered phase which is neither compressible nor Mott 
insulating. This state is characteristic of quasidisorder in the presence of a strong trapping potential. 
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I. INTRODUCTION 

The study of quantum phases for the strongly cor- 
related system of identical particles in a lattice is a 
core research topic in condensed matter and solid state 
physics Due to rapid advances in the field of ul- 
tracold atomic physics such a system has been routinely 
realized in optical lattice experiments, mimicking con- 
densed matter physics with unprecedented precision 
Indeed, ultracold atoms in optical lattices are considered 
as promising quantum simulator for the strongly corre- 
lated electron system in solid state physics Some of 
the hallmarks of such trapped quantum gases in a lattice 
are the experimental realization of the supcrfluid-Mott 
insulator transition Q, of a TonksGirardeau gas [H, of 
the Anderson localization for a Bose-Einstein condensate 
with tunable interactions [(|, and many others. 

In the experiments ultracold systems are usually real- 
ized in the presence of external parabolic confinement, 
hereafter, we synonymously call it trapping potential. 
This makes the ultracold atomic system spatially local 
since the local chemical potential varies at each spatial 
point due to inhomogeneity. The competition between 
the atom-atom interaction energy and the local confin- 
ing potential at each point dictates the true state of the 
system. In particular, in a strongly correlated system 
with random or quasiperiodic disorder the trapping po- 
tential may have a dramatic effect in the ground state 
properties and thus demands a detail and systematic in- 
vestigation. In fact, there have been a number of studies 
that show the importance of the inhomogeneity in the 
lattice system 0. 

A Mott state for the trapped fermionic atoms in an 
optical lattice and its incompressible nature has recently 
been studied experimentally @, The system's negligi- 
ble response with the variation of the trapping potential 
is the characteristic feature of the incompressible Mott 



insulating state. A similar line of thought has also been 
proposed in the context of detecting Bose-glass phase in a 
lattice system with random or quasi-random lattice po- 
tential [l(j-[l2|- This technique the so-called the trap- 
squeezing spectroscopy (TSS) relies on the response of 
the system with the variation of the trap pin g potential, 
since the size [l(| and the central density jl2| of the sys- 
tem depend on the strength of the trapping potential. 

In Ref. [l(| it is reported the fact that the variation 
of the radius of atomic cloud with the confining poten- 
tial may be a useful tool for detecting the disordered 
phase such as the Bose-glass [l3| . It should be noted that 
for a trapped system there is a nonvanishing supcrfluid 
fraction since the atoms near the edges can be delocal- 
ized with significant number fluctuations. In this sense 
it is more precise to call the superfiuid (Mott) phase in 
a trapped system as the phase with high (low) degree 
of compressibility. Finding a true signature of incom- 
pressible region over many lattice sites requires a strong 
lattice with minimal edge fluctuations. The Mott and 
the superfiuid region can be dramatically affected by the 
variation of the external trapping potential as well as the 
total number of atoms present in the system since the 
chemical potential depends on both of them. Fixing the 
chemical potential by varying the atom number when the 
trap frequency is changed has also been proposed in fTlj j . 
However, in real experimental situation this restriction 
poses two challenges: (i) the chemical potential is not 
directly accessible in an experiment, and (ii) the number 
of atoms in each experimental shot is fixed. 

In this paper we follow a more direct approach. 
We study the properties of strongly interacting bosonic 
atoms in an optical lattice in the presence of external 
trapping and additional quasiperiodic potential. For a 
strong lattice the tunneling is exponentially suppressed so 
that the system can be in deep Mott regime with vanish- 
ing number fluctuation. We first develop a phase diagram 
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in the parameter space spanned by the total number of 
atoms and the external trapping potential. In a pure 
system the phase diagram consists of alternate stripes of 
incompressible Mott and compressible superfluid region. 
Such compressible-incompressible region can also be ob- 
tained by studying the variation of the root mean square 
(RMS) size of the cloud with the number of atoms as well 
as with the trapping potential. In the Mott state the size 
of the cloud varies linearly with the atom number since 
the available state for additional particle is located at 
the edge of the cloud. For a small number of atoms the 
incompressible Mott state is robust with respect to the 
variation of the trapping potential. This state becomes 
compressible when the disorder is introduced by means of 
an additional lattice potential. In the quasiperiodic sys- 
tem with the fixed number of atoms, the decrease in cloud 
size with increase in trapping potential along with the 
vanishing superfluid fraction suggests that the system is 
in compressible Bose-glass phase. In addition to the com- 
pressible state, the RMS size of the cloud also develops 
plateaus as the trapping potential becomes stronger indi- 
cating the presence of incompressibility in the disordered 
system. This state differs from the Mott phase particu- 
larly in their density distribution, and is characteristics 
of the quasiperiodic potential in the presence of trapping 
potential. Similar findings but in different model in the 
presence of off-diagonal disorder has also been reported 
using quantum Monte-Carlo simulations (l4| . 

The paper is organized as follows. In Sec. II we intro- 
duce the Bose-Hubbard model in a one-dimensional lat- 
tice in the presence of additional modulating potential. 
We employ decoupling mean-field method (DMF) [l5[ 
to solve the ground state of the Hamiltonian. For com- 
pleteness, we also connect the Bose-Hubbard parameters 
approximately to the experimentally measurable quanti- 
ties such as lattice heights, atomic scattering lengths etc. 
In Sec. Ill we study the properties of the system partic- 
ularly the quantum phases as the trapping potential and 
the atom number are varied. Finally in Sec. IV we study 
the effect of quasiperiodicity, in particular, the possible 
modifications in quantum phases and the appearance of 
disordered Mott-likc incompressible phase. 



II. QUASI DISORDER BOSE-HUBBARD 
MODEL 



At zero temperature, the quantum state of bosonic 
atoms in a one-dimensional deep optical lattice with diag- 
onal disorder is approximately governed by the modified 
Bose-Hubbard (BH) Hamiltonian Q3, [HI 



- J(aja i+ i +h.c.) + — hi(hi-l) - ^hi + eitii 



(1) 

where d\ and di are the usual creation and annihilation 
operators for a boson at the i th lattice site, and hi = djd; 



is the corresponding number operator. Here, J is the 
kinetic energy for an atom hopping to the nearest site, 
U is the on-site two-body interaction energy, fi is the 
global chemical potential, and is the on-site energy 
of an atom. The latter is composed of both harmonic 
potential resulting from magnetic and optical fields and 
disordered potential generated by a secondary lattice [13, 
El, 



VtZ 2 , A _2 



A cos (irai + S) . 



(2) 



Here the trapping strength is Vt = mu> z d , m being 
the mass of an atom, uj z the longitudinal trapping fre- 
quency and d the lattice spacing. The incommensurate 
parameter a is the ratio between the wavelengths of the 
primary and secondary laser lights, whereas A, which 
is proportional to the intensity of the secondary lattice, 
measures the strength of the disorder . The relative 
phase between the primary and the secondary lattices 8 
is chosen in such a way that the harmonic potential min- 
imum lies on one of the minima of the secondary lattice 
potential, preferably chosen at the central lattice site. All 
the parameters discussed in this model can be expressed 
in terms of experimentally measurable quantities such as 
the heights of the primary and secondary lattices, the 
wavelengths of the lasers, the s-wave scattering length 
characterizing the atom-atom interactions etc. (l9l . |20| 
(see Sec. III). 

In a homogeneous case the Hamiltonian ([T|) shows the 
superfluid-Mott insulator transition when the ratio J/U 
exceeds some critical value (J/U) c 0, [U|. In a weakly 
interacting regime, the tunneling term in the Hamilto- 
nian dominates over interactions, and the corresponding 
many-particle ground state is the superposition of all sin- 
gle particle wave functions spread over the entire lattice, 
causing the system to be superfluid. On the other hand, 
if the tunneling term is negligible the ground state en- 
ergy of the interacting system is minimized with a well 
defined number of atoms per site and the many-body 
ground state is a product of Fock states at each lattice 
site leading to the Mott insulating state. In this regime 
the spectrum presents a gap of the order of U, corre- 
sponding to the displacement of an atom to an occupied 
site. 

In addition to these phases the BH model also predicts 
an insulating but gapless phase in the presence of random 
disorder, known as Bose-glass However, the search 
and detection of this phase has been a long standing ex- 
perimental challenge from the time of its conception, and 
so far there are preliminary experimental indications (l8j 
but no clear evidence that directly supports its existence. 

The quantum properties of lattice bosons governed by 
the BH model can be studied accurately through exact 
numerical diagonalization only for a few lattice sites and 
a small number of particles [21(. Other methods, such 
as quantum Monte Carlo (QMC) simulations Ell |22| . 
density matrix renormalization group (DMRG) [23| and 
time-evolving block decimation (TEBD) [24[ have been 
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extensively used to extract the information on the ground 
state and excitation properties of the system. In this pa- 
per we use the decoupling mean-field (DMF) approxima- 
tion [l5j , that provides a reasonable picture of the phase 
diagram of the homogeneous system for both pure and 
disordered case. This method also parallels to another 
popular mean-field approach the so-called the Gutzwiller 



approximation [251 ]. 

In the DMF approach one can approximate the off- 
diagonal terms in the Hamiltonian as 



(3) 



yielding the mean-field Hamiltonian 



J 



Hmf = ^ 



J((f>*a i+ i + a}(f>i + i - 4>*<t>i+i + h.c.) + ^-hi(hi - 1) - /m, + e,n 



(4) 



r 



where <pi = (di) is the superfluid order parameter, that tonian can be written as the sum of the on-site hamilto- 
provides a residual coupling between neighboring lattice mans as Hmf = J2i Hi with 
sites. By using periodic boundary conditions the Hamil- 



Hi = J{4>t- X + 4>t+i)ai + J{4>i-i + 4>i+i)a\ - J(0i_i + , 

I 



■> i+ i)4> t + ^hi{hi - 1) - [ihi + eihi 



(5) 



The Hamiltonian can be diagonalizcd self-consistcntly us- 
ing the number basis and the ground state can be found 
for a given set of parameters J, U : and A [l5| . 



A. System geometry and parameters 

In recent experiments it has been considered a system 
of ultracold bosons in the presence of two strong trans- 
verse optical lattices of intensity sj_ in units of photon 
recoil energy, creating a two dimensional array of inde- 
pendent tubes Q. For large values of s±, the coupling 
between different tubes vanishes, so that each tube por- 
trays a quasi ID atomic system. 

In order to study the feature governed by the Hamilto- 
nian (JTJ a lattice can be created in the tube direction [z- 
axis) by shining a third pair of counter-propagating laser 
of wavelength Ai of lower intensity, s\ <C s±. Quasidis- 
order can be incorporated by introducing an additional 
pair laser of much smaller intensity sf and of wavelength 
A2 incommensurate with the main lattice, producing an 
effective quasiperiodic potential [l8[. Such a controlled 
disorder essentially captures the features of random dis- 
order as demonstrated in a recent experiment on Ander- 
son localization [fjj], and is therefore a suitable medium 
to study disorder induced quantum phases such as Bosc- 
glass. 

Under these conditions, the BH parameters can be 
tuned with the experimental measurable quantities, ac- 



cording to [19|, l2f 
J 
U 
A 



1.43 



( S D a98 exp(-2.07 v ^f), 



5.97(a s /Ai)( S f 4) a88/3 . 

fM! P -/37W) 7 



(0) 
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(8) 



in terms of the recoil energy En = 2h 2 ir 2 /m\ 2 . Here a s 
is the s-wave scattering length and a, f3 and 7 are fitting 
parameters. Each laser beam also creates an additional 
harmonic confinement (due to the focusing) that intro- 
duces a source of inhomogencity; therefore the number of 
atoms in each tube varies with maximum density at the 
central tube. In the experiments these numbers also de- 
pend upon the loading procedure, the adiabaticity of the 
lattice ramping and the external harmonic confinement. 
If the cloud follows the Thomas-Fermi (TF) profile, one 
can write the atom distribution in each tube as 1271 



N~ v = N 0fi 1 - 



2\ 3/2 



B 2 



(9) 



where (cc, y) represents the position of the tube in the 
transverse plane, Rtf is the TF radius given by 



Rtf 



tvjj 



d 2 



15Na s d 2 

a-ho 



dxf (x) 



2-1 1/5 



, (10) 



with N being the total number of atoms, and iVo.o = 
hN / (2ttR^ f ) the number of atoms in the central tube. 
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Here m, a s and d are respectively the mass, the s-wave 
scattering length and the lattice spacing in the direction 
of the tube, whereas wj_ and Q = (ui 2 1 _u) z ) 1 f i are the trans- 
verse and averaged trapping frequencies, a\ m = \J h/muj z 
is the harmonic oscillator length, and fo is a function 
localized at the origin. 

For the discussion presented in this paper we choose 
the typical parameters of the experiment in Ref. fl8| . 
In particular, for a system of N — 2 x 10 5 87 Rb atoms 
with s-wave scattering length a s = 100 ao, in trapping 
frequencies ui± = 2irx 90 Hz and uj z = 27rx 77 Hz, and the 
transverse lattice height s± = 40, the number of atoms 
in the central tube and the average number of atoms in 
each tube are respectively -/V .o ~ 200 and N av rj 80. 

Since our model is strictly one-dimensional we fix the 
relevant parameters as stated in the previous paragraph 
by paying attention to a single tube only. In all our 
calculations the height of the primary lattice is fixed to 
sf = 16 corresponding to the ratio U/J ~ 132. The 
length of the lattice is fixed to L = 64. We study the 
properties of the system by varying both the number of 
atoms in a given tube and the trapping potential. 



characterized by considering several quantities, including 
condensate fraction, supcrfluid fraction, number fluctua- 
tions, defined as follows. 

The condensate fraction is defined as the ratio between 
the number of condensate atoms to the total number of 
atoms, and can be simply expressed in terms of the order 
parameter <f>i as 



fc 



Ef Ifr 

Aft 



(11) 



where Aft is the number of atoms in the tube. The su- 
perfluid fraction, on the other hand, can be calculated by 
measuring the response of the system to an infinitesimal 
flow, by imposing twisted boundary conditions [U [29[ 
that corresponds to a phase gradient across each lattice 
site. The supcrfluid component responds to the velocity 
field without excitations, and the gain in the kinetic en- 
ergy is mainly due to the superflow of the condensate [28| . 
The superfluid fraction can be written as [29| 



2mL 2 Et, 
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III. PRIMARY LATTICE ONLY: PURE CASE 

In this section we analyze the behavior of the system 
with the variation of the confining potential and the to- 
tal number of atoms for the case when there is no addi- 
tional secondary lattice. The state of the system can be 



where 9 is the the phase twist along the lattice and Eg , Eq 
are the resulting ground state energies with and without 
the twist. This difference in energies with and without 
phase twist can be treated perturbatively for small 6. 
Within the second order perturbation theory we get the 
expression for the superfluid fraction as 



iEi 
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+ 4>i+i)4>i 



j L 



-<^+i) : 



En V^ + T/:+l/n°- 



(13) 



From this expression it can be seen that the compe- 
tition between the first and the second order terms give 
rise to the net supcrfluid fraction, the role of the second 
order contribution is to damp and reduce the supcrfluid 
component. It is easy to verify that in the homogeneous 
case, the condensate and superfluid fractions coincide, 

fs fc- 

Atom number fluctuations at each lattice site are de- 
fined as 

An, = ^(n 2 ) - (h,) 2 , (14) 

and directly characterize the different quantum phases 
of the system, e.g., in the Mott phase the number fluc- 
tuation vanishes Ani/ni = 0, whereas in the supcrfluid 
phase it becomes large, An, ~ nj. 

In Fig. [T]wc show a typical ground state density dis- 
tribution for a trapping strength V t = 0.001 (that cor- 
responds to the longitudinal trapping frequency ui z sa 



2n x 70Hz) and different values of the total number of 
atoms, Aft, in the tube. The atoms are delocalizcd with 
significant supcrfluid fraction for a few number of atoms 
in the tube. As the number of atoms reaches the critical 
value a flat distribution with Mott-like structure appears 
near the center. Since the Mott phase is incompressible, 
further addition of particles uniformly expand the cloud 
until a non-integer occupancy at the center (n c ) = n + e 
is energetically favorable. This gives rise to large number 
fluctuations and supcrfluid fraction at the center of the 
trap. 

Figure [5] depicts the density plot of the superfluid frac- 
tion in the Aft - Vt plane for the parameters discussed 
in Sec. IIA. The brighter region represents the nonva- 
nishing superfluid fraction while the darker region cor- 
responds to the insulating state. Each successive su- 
pcrfluid stripe represents the state of the system char- 
acterized by the density at the center of the trap with 
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FIG. 1: (Color online) Atom distribution in the lattice for the 
trapping strength Vt/Eu = 0.001 as the number of atoms is 
increased from Mt — 2 to Mt = 74. The length of the lattice 
is L = 64 and the height of the primary lattice is sf = 16 
corresponding to the ratio U/J = 132. 



FIG. 3: (Color online) Cloud size as a function of the total 
number of atoms for different trapping potentials Vt in units 
of Er. The Mott phase is incompressible as the size of the 
cloud for different trapping potentials for a fixed atom number 
remains constant. The wide separation between the curves 
represents the compressible superfluid phase. 
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FIG. 2: (Color online) Density plot of the superfluid fraction 
in the Nt - Vt plane for the same parameters of Fig. [T] Each 
stripe represents the nonvanishing superfluid fraction. 



(n c ) = 1 + e, 2 + e, 3 + e etc. These stripes bend to- 
wards the smaller number of atoms when the strength of 
trapping potential is increased. 



A. Size of the cloud and compressibility 

A direct manifestation of the incompressible (com- 
pressible) nature of the Mott (superfluid) phase can be 
obtained by measuring the response of the cloud as the 
trapping frequency is varied. The system is incompress- 
ible if the size of the cloud for a fixed atom number re- 
mains unchanged. This method has been successfully 
implemented experimentally to detect the Mott phase in 
a fermionic system || Q . 

The global compressibility can be related to the size of 



/, . , dlZ 



(15) 



with 1Z being defined as the root-mean-squared (RMS) 
width 



\ 



(16) 



in units of lattice constant d. The latter can be measured 
experimentally with great accuracy by means of in situ 
imaging technique [31 1. 

In Fig. [3]we show the RMS size 1Z of the atomic cloud 
as a function of the total number of atoms Nt for dif- 
ferent trapping potentials. In the Mott regime, the size 
of the cloud increases linearly with the number of atoms 
until the potential energy of an atom at the periphery 
exceeds the onsite interaction energy A sharp kink in 
the value of 1Z in each curve represents a transition from 
occupancy (hi) = n to (hi) = n + e at the center of the 
trap. This transition occurs earlier for the higher value 
of the trapping potential. 

Once the system exceeds the critical occupancy the size 
decreases with the addition of atoms since more nearby 
central sites allow multiple occupancies, and therefore 
the RMS size decreases according to Eq. (fT6|) . A wide 
separation between these curves for close values of Vt 
indicates that the system is compressible. In the Mott 
regime, these curves overlap. This is the case of the Mott 
phase for low filling in Fig.[T] corresponding to the initial 
linear slope indicated in Fig. [3J 

This can also be seen in Fig. 2] where we plot the RMS 
size as a function of the trapping potential Vt for different 
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FIG. 4: (Color online) Cloud size as a function of the trapping 
potential for different atom numbers Nt = 25, 40, 60 and 80. 



atom numbers Nt = 25, 40, 60 and 80. For a small atom 
number the size virtually remains unaffected, whereas it 
decreases for larger number of atoms. 

The slope of each curve is proportional to the com- 
pressibility (sec Eq. (|15[1 ). The system is incompressible 
for low atom number and for low trapping potentials, 
when the atom distribution is independent on the trap 
strength. This situation is the boson counterpart of the 
fermionic Mott insulator as demonstrated in [1, Q . 

For larger number of atoms the distribution spreads 
over the outer part of the trap and therefore, an increase 
of the trapping potential produces a large energy offset 
with respect to the center of the trap and atoms jump 
towards the center of the trap resulting a dramatic drop 
in the cloud size. A similar observation has also been 
pointed out in Ref. using TEBD calculation. For 
larger number of atoms far from Mott filling, multiple 
occupancy smoothly fill the energy gaps and the system 
always responds to the trap variation resulting a com- 
pressible superfhiid phase. 

The mechanism of supefluid-Mott transition by vary- 
ing either the number of atoms or the harmonic confine- 
ment can also be explained in the phase diagram shown 
in Fig. [2] One can move along the vertical direction for 
a fixed number of atoms or along the horizontal direc- 
tion for a fixed trap strength to see if there is a dramatic 
increase in the supcrfluid fraction. Each bump in the 
supcrfluid fractions represents the jump in atom number 
from an integer to non integer value in the center of the 
trap. 



B. Plateaus and emergence of mini gaps 

In Fig. [5] we show the behavior of the total number 
of atoms as a function of the chemical potential for dif- 
ferent external trapping potentials for the same param- 
eters considered before. In the Mott regime the curve 
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FIG. 5: (Color online) Variation of the total atom number 
Nt as a function of chemical potential /i for different trapping 
potentials Vt ranging from Vt — to Vt — 0.003 in units of 
En. The Mott plateaus gradually shrinks as Vt gets stronger 
and finally disappears when Vt > Vt - For strong confinement 
the curve develops a series of plateaus (inset) which becomes 
more pronounced when fi/U — > 1. 



shows plateaus for integral fillings. In addition to the 
Mott plateaus the interplay between lattice periodicity 
and trap inhomogeneity displays a series of additional 
plateaus in the density versus chemical potential curve 
(see inset). Such plateaus become more pronounced when 
the trapping potential is strong. The origin of these 
discrete jumps in the density is due to the energy off- 
set produced by the harmonic potential at each lattice 
point. For strong atom-atom interactions the system 
avoids multiple occupancies and the corresponding en- 
ergy is minimized by filling the unoccupied lattice site 
starting successively from the center of the trap analo- 
gous to fermions filling in the Fermi sea. As the number 
of atoms increases the addition of an extra atom to the 
system already contained Nt atoms costs energy of the 
order of j|mw 2 (A/t/2) 2 <i 2 which is the energy gap as seen 
in Fig. [5] (inset). As expected this gap becomes more 
pronounced when the trapping frequency is strong. 

In addition to these plateaus Fig. [5] also shows the dis- 
appearance of the Mott plateau as the trapping potential 
is reached to a critical value V£. For this particular set of 
parameters the critical trapping potential is V t c f=s 0.001. 
For smaller trap strength, the strong interaction pushes 
the atoms to the boundary so that the system reaches 
integer filled homogeneous Mott state. As the trapping 
potential is increased the true Mott region with (n) = 1 
shrinks until double occupancy is more favorable and the 
atoms start jumping to the central lattice site. Unlike in 
the homogeneous system, the nonzero average slope in 
the N versus \x plot for each value of the trapping po- 
tential gives a small but finite compressibility and hence 
the trapped system is, in principle, always compressible. 
Also note that the compressibility k tends to vanish when 
fj,/U 1 for each value of Vt > V t c since k — dn/dfj,. 
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FIG. 6: (Color online) Cloud size as a function of the total 
number of atoms for different trapping potentials for the dis- 
order strength A = 0.5 and A = 0.01 (ss 2 J) (shown in the 
inset). Both A and Vt are expressed in terms of the recoil 
energy E R . 



FIG. 7: (Color online) Cloud size as a function of the trapping 
potential for the total atom number in the tube Aft = 25. Also 
shown (in the insets) are the atom densities in the weak (a) 
and strong (b) confinements. 



IV. EFFECTS OF SECONDARY LATTICE: 
DISORDERED CASE 

In this section we consider the situation when an ad- 
ditional secondary lattice, incommensurate to the main 
one is introduced resulting a quasiperiodic system as dis- 
cussed in Sec. II. Quasiperiodic disorder can be con- 
trolled by varying the intensity of the secondary laser 
as the strength, A, is proportional to the height of the 
secondary lattice sf (see Eq. flSJ). The lattice incom- 
mensurate parameter gives an additional characteristic 
periodicity of 1/(1 — a). 

In the following we will consider a = (v5 — l)/2, the 
inverse of the golden mean [30| . In a non- interacting sys- 
tem the threshold for the extended to the localized single 
particle state occurs at A/ J = 2 [ljj. In the strongly 
interacting regime, J <C U, the interaction is the only 
dominant energy scale. Therefore, the localization may 
occur as soon as the disorder is switched on [13| . 



A. Size of the cloud and compressibility 

The effect of trapping potential on the size of the cloud 
in the presence of the secondary lattice is demonstrated 
in Fig. |B1 where we plot the RMS size 71 as a function 
of the total atom number Aft , for two values of the disor- 
der strength namely A = 0.5 (3> J), which is the critical 
value of the disorder strength for the disappearance of 
the Mott lobe EH, and A = 0.01(w 2 J) for which the 
localization-delocalization transition in the noninteract- 
ing system occurs [T3|. The rest of the parameters are 
the same as in the pure case. The robust nature of the 
Mott insulating state for small A = 0.01 (shown in in- 
set) and a wide separation in TZ for large A = 0.5, as 



the trapping strength is varied, manifestly indicate that 
the former state is incompressible while the latter state 
is compressible. It is noted that there are certain points 
for A = 0.5 at which the curves for different values of 
the trapping potentials tend to cling together when the 
number of atoms are small. This behavior may be ac- 
counted due to the presence of additional plateaus as 
seen in Fig. [SJ 

In Fig. [7] we plot the cloud size TZ as a function of 
trapping potential for a fixed atom number Aft = 25, 
corresponding to the Mott state in the pure system. For 
small trapping potential the size TZ decreases monoton- 
ically. We also show the typical state of the system for 
the two values of the trapping potentials (in the insets). 
The atom density is inhomogeneously distributed due to 
the quasipcriodicity with many unoccupied lattice sites 
in between. As the trapping potential increases more and 
more atoms from the edges displace towards the center 
of the trap or near the quasiperiodic minima and thus 
reducing the size of the cloud. This decrease in TZ along 
with the vanishing superfluid fraction may hint that the 
system has entered into the compressible (quasi) Bosc- 
glass phase. This prescription works fine as long as the 
trapping potential is weak. 

As the trapping potential gets stronger, atoms start 
accumulating near the center of the trap and the size 
TZ develops intervening plateaus, indicating incompress- 
ible nature of the system. These plateaus becomes more 
pronounced when the trapping potential is strong. From 
the inset (b) of Fig. [7]we see that the state of the system 
does not change over a wide range of trapping poten- 
tials. This incompressible state in the presence of disor- 
der differs from the Mott phase particularly in its density 
distribution. The atom distributions are still inhomoge- 
neous with many intervening unoccupied along with the 
doubly occupied sites. These insulating states are the 
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characteristics of the quasiperiodic potential in the pres- 
ence of trapping potential and resemble to the so-called 
Mott-glass phase as reported in Ref. fl4j |. 

V. CONCLUDING REMARKS 

In conclusion we have studied the possibility of detect- 
ing quantum phases of strongly interacting bosons in a 
one dimensional optical lattice in the presence of trapping 
and quasiperiodic potentials. In general the trapping po- 
tential may obscure the superfluid-Mott transition and a 
direct mapping of the phase diagram from the inhomo- 
geneous to the homogeneous system can be challenging. 
In particular implementing the idea of trap-squeezing to 
detect the quantum phase poses challenges due to the 
local nature of the system [12| ■ 

In this paper we have studied the effects of trapping po- 
tential in a strongly interacting bosonic system by vary- 
ing the number of atoms. In a pure system a true incom- 
pressible Mott state can be observed for small number of 
atoms over a range of trapping potential. In parallel to 
the recent experiments on detecting the incompressibility 



of fermions [1, Q , the robust nature of the Mott state is 
verified by varying the trapping potentials as the cloud 
size remains unaffected. However, the introduction of 
the disorder induces trap-dependent cloud size. The van- 
ishing supcrfluid fraction along with the trap-dependent 
cloud size may hint that the system is in Bose-glass phase. 

In addition to the compressible phase, the strong trap- 
ping confinement may induce a Mott-like incompressible 
states in the presence of disorder. These insulating states 
differ from Mott phase in their density distribution. The 
strongly interacting system in the presence of disorder 
interpolates between two insulating states namely incom- 
pressible Mott-like states and compressible Bose-glass as 
the trapping potential is continuously varied. 
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